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Abstract 


A model is presented to enplain the highly variable yet low level of 
Langmuir waves measured in situ by spacecraft when electron beams associated 
with Type III solar bursts are passing byi the low level of excited waves 
allows the propagation of such streams from the sun to well past 1 AU without 
catastrophic energy losses. The model is based, first, on the existence of 
large scale density fluctuations that are able to efficiently diffuse small-k 
beam-unstable Langmuir waves in phase space, and, second, on the presence of a 
significant isotropic non-thermal tail in the distribution function of the 
background electron population, which is capable of stabilizing larger k 
modes. The strength of the model lies in its ability to predict various 
levels cf Langmuir waves depending on the parameters. This feature is 
consistent with the high variability actually observed in the measurements. 

The calculations indicate that, for realistic parameters, the roost unstable, 
small k modes are fully stabilized while some oblique mode with higher k and 
lower growth rate might remain unstable j thus a very broad range of levels of 
Langmuir waves is possible from levels of the order of enhanced spontaneous 
emission to the threshold level for nonlinear processes. On the other hand, 
from in situ measurements of the density fluctuations spectrum by ISEE 1 and 2 
in the viclrity of the earth, it is shown that measured 100 km scale 
fluctuations may be too effective in quenching the instability. If such 
strong density fluctuations are common in the solar wind, we show they must be 
highly anisotropic in order to allow the build-up of Langmuir waves to the 
observed mV/m range. Moreover, the anisotropy must be such that the strongest 
variations of density occur in a plane perpendicular to the magnetic field. 


1. Introduction 


Highly variable yet rather low levels of Langmuir waves have been 

I 

measured in situ by spacecraft when electron-beams associated with Type III 
solar bursts are passing by, due to the low level of excited waves the 
propagation of such streams from the sun to well past 1 AU is possible without 
catastrophic energy losses. Detailed measurements at 1 AU (Lin et al., 19B1) 
have shown clearly the simultaneous occurrences of a bump on tail electron 
distribution function and a rise of plasma waves above the background level 
Yet they also reveal that this distribution remains for a long enough time 
(t 10 min) to drive the plasma waves to far higher levels (according to linear 
instability theory) than the maximum values observed of a few mV m"^ and 
that, contrary to what might be expected from pure quasilinear theory, no 
plateau formation of the distribution is evident. One possible explanation is 
that various nonlinear wave-wave interactions such as decay instability, OTS. 
modulational instability, and soliton collapse are effective in limiting the 
wave growth by shifting the waves out of resonance with the beam (Ueatherall 
et al., 19Q1; Grognard, 1962). It should be noted, nevertheless, that the 
threshold for such nonlinear processes to occur is about 10 mV m“^ for 
conditions prevailing at 1 AU. 

On the other hand, the electron density has been observed to fluctuate 
significantly in the solar wind. In situ measurements of the density 
fluctuation spectrum by ISEE 1 and 2 (Celnikier et al., 1983) reveal that a 
variation 6n/n as high as 10 ^ may exist on the 100 Km scale range in the 
vicinity of the Earth, while observations of interplanetary scintillation from 
extragalactic radio sources (Cronyn, 1972) lead to an average value 6n/n of 


P G 1 


the order of 10"^ (Smith and Sirae, 1979). Thus it seems that the 
background density of the solar wind, although not precisely known, is far 
from being homogeneous and that, when dealing with the Langmuir turbulence 
associated with Type III bursts, typically at Km wavelength, the effect of the 
inhomogeneities must be accounted for. In this paper we show that the large 
scale density fluctuations may be extremely effective in shifting the waves 
out of resonance with the beam and therefore in quenching the instability. 
Following an idea studied by Nishikawa and Ryutov in the context of laboratory 
relativistic electron beams (Nishikawa and Ryutov, 1976) and by Goldman and 
DuBois in that of the solar wind (Goldman and DuBois, 1982), it is 
hypothesized that the effect of the inhomogeneities on the Langmuir turbulence 
may be treated in terms of a diffusion process in phase space. However, 
unlike the latter authors, we take into account the non-Maxwellian background 
electron distribution, which results in quite different spectra, and we do not 
assume that the background density fluctuations are necessarily isotropic. In 
fact, it will be shown that, if the high level of density fluctuations 
measured by ISEE is common and of long duration, the fluctuations must indeed 
be highly anisotropic in order to allow the build-up of Langmuir waves up to 
the mV/m range. Moreover, the anisotropy must be such that the strongest 
variations of density occur in a plane perpendicular to the magnetic field. 
This feature leads to a picture of the solar wind as composed of structures 
elongated along the magnetic field axis in a way consistent with the ducting 
considerations invoked by experts in radio-wave propagation (Robinson, 1983 j 
Duncan, 1979s Bougeret and Steinberg, 1977). 


2. Model 


A heuristic picture of the wave diffusion process nay be outlined briefly 
for the reader familiar with particle velocity space diffusion (i.e. quasilinear 
Vlasov theory). In the Langmuir wave diffusion, the plasmon (of wavenumber k) 
plays the role of a diffusing particle in the usual quasilinear Vlasov 
theory and the role of the external potential is played by the varying 
refractive index associated with the density fluctuations (with a spectrum of 
wavenumbers characterized by q). Suppose q~^ < L^ where L^ is the 
typical growth length of the Langmuir waves so that a spectral energy density 
averaged in space U(k) may be used. Further, we assume q << k; then the 
relative change of k is small in each scattering event and the process is 
treated as a diffusion. The diffusion in k-space experienced by the plasmon 
may be estimated as D k t^^ with ^ the fluctuating force 

acting on it for an autocorrelation time t^^. Since ^ x sn/n q, 

and tg^ is of the order of the transit time for the plasmon across one density 
fluctuation, t^^ s CqVg)“^ s (qkXjjVg)’^ , we have D = (q/k) (kx^)’^ (Sn/n)^. 

A more rigorous general theoretical derivation (Goldman and DuBois, 1982) 
leads, with the above assumptions, to the following Pokker-Planck equation for 
the spectral energy density of Langmuir waves in phase space. 

®t ”k “ 2^k“k ^ ®kj °j£ ®k^ ”k * h 

where the diffusion coefficient is 

= «pe " J C(a) q^q^ 8(2k-a) 


( 2 ) 


and C(£) is the spectrum of density fluctuations i 


/ d^q/C2Ti)^ C(fl) « (Bn^/n^) . 


( 3 ) 


The diffusion coefficient, Equation (2), incorporates the assumption of elastic 
scattering. That assumption is justified when the propagation velocity of the 
density inhomogeneities is small as compared to the group velocity of the plasmon 
Vg = 3 l‘^pVg. Note that in Equation (1) we have added an electron spontaneous 
emission term S|^. 


By assuming that the problem displays an azimuthal symmetry around the 
beam axis we may rewrite Equation (1) as 
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3j. W(k,e^ = 2 yjj W{k,0) + (sine) '’8^ D(k,6) sine 8gW(k,0) + S(k,0) 


(4) 


where D(k,e) is derived from Equation (2): 

Dfk,0) «= 1/k^ 1/(12 ti^) fl(0) 


(5) 


with 

cos^e/sin^e 

Qfe) = J dCcoB 0) c(§) 5 r;;; r 

(sin'^esin^e-cos^ecos 


( 6 ) 


where C(e) = J* dq q^ C(q,0) and I is the domain cose <|sin0l. 

Note that the diffusion occurs in angle only as a result of the assumption of 
elastic scattering. 
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The growth rate (or damping) and the spontaneous emission term Sj, 
in Equation ( 4 -) depend on the actual distribution function. We model them by 

f(v) = n/f2w exp (- v^/2v„^) 

C B * 

+ 5v^^/4ti v'® HCv-v^) 

+ ov‘^ exp [- (y-Vj^l^/2ov^] . 


Here H is the Heaviside step function, n,p and nj^ are the fractions of 
electrons in an observed non-thermal electron tail and in the beam feature 
(relative to the bulk); 6v is the beam velocity spread. The shape of the 
tail simulates a non-thermal component, such as that observed by Lin et al. 
(1901) prior to and after the passage of the beam. We assume that it is valid 
down to a velocity v - v,p b 5v^. Since the beam velocity v^j s 60-90 Vg, 
waves with a phase velocity v^j^ < 5Vg will not be considered. Then the 
presence of a step at in the distribution function has no importance, and 
a smooth matching of the tail to the bulk is not required. It should be 
mentioned here that attempts have been made to simulate the slight measured 
anisotropy of the tail associated with a heat flux (Gurnett et al., 1979); 
these have not led to significantly different results of concern to us, so 
that complication will be omitted in the following. 

The growth/damping rate consists thus of two isotropic terms due to 
the Landau damping from the bulk velocity distribution (denoted by M for 
Maxwellian) and from the tail (denoted by T) 


= _ (Ti/0)^''^ (lAXp)^] 


(0a) 
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^ ‘ ® “ (Bb) 

and of an anisotropic term due to the beam 

B cose-w /k)exp[-l/2Av^(v. cose-« /k1^] 

V V o pe 

Similarly, the spontaneous emission is composed of tno isotropic and one 
anisotropic term 


sfM) = mv2 . 


f9a^ 


= (5n/121 w"g mv^ n^. (kv^l^ , 


(9b) 


3^*=^ B (tt/ 2 ) 1/2 „2^ mfWpg/k)2 n^/kav exp l/26v^ (Vj^cose - Wpg/k)^] 


(9c) 


When the rate of diffusion of the Langmuir waves out of unstable regions 
into stable regions of the phase space is high enough, Equation (4) has a 
stationary solution. In order to study the stability boundaries in parameter 
space as well as the level of electric field, Equation (4) has been solved 
numerically in the steady state limit for various parameters. 
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3. Results and Interpretation 

3.1 Numerical Results 

The Fokker-Planck Equation (4) displays two kinds of solution. Either 
the rate of diffusion of the Langmuir waves out of unstable regions into 
dissipation regions of the phase space is rather small and the wave energy 
grows up ad infinitum as in the usual linear theory, or the rate is high 
enough in order that the wave energy saturates, and Equation (41 then has a 
steady solution. 

A typical steady solution may be seen in Figure 1. Figure la displays the 
contours of the wave energy density spectrum Wj^/Tg in the k-e plane of 
phas^ space. Note that the enhanced spontaneous emission due to the 
non-thermal tail, Wj^/Tg a , would lead (with the 

parameters of the Figure) to a spectrum of order 10^ in the region of phase 
space displayed. Thus, because of the diffusion, even the backward 
propagating damped modes (0>n/2) have a comparatively high level of wave 
energy. Let us now consider the total electric field assoc;, ated to a 
steady solution. 

3.2 Properties of Steady State Solutions 

On a sphere of radius k in three-dimensional k-space, energy source 
regions are connected to energy sink regions by angular diffusion. Therefore, 
it is possible to write a balance equation after integrating Equation (4) on the 
sphere. For that purpose we define the two averages 

- 

<Yj^> = 1/Wjj J ^ d(coE0) Y|^(0) W(k,0) 


(101 



where 


« J ^ d(cos0) W(k,e) 


and 


= 1/2 J, d(coB0l 

« -1 (in 

Since the angle average of the diffusion term vanishes, the balance equation 
may be written as 




( 12 ) 


Where the diffusion coefficient D(k.e) is large, a very weak angular 
gradient in the spectrum is sufficient for the diffusion fluK to remove the 
wave energy from regions of sources into regions of sinks. The spectrum is 
therefore approKiraately Independent of e so that <y^> = Vj^. A key 
point of Nishikawa and Ryutov (1976) was to show that for any 
distribution irrespective of the presence of a beam. Where D(k,e) is weaker, 
the spectrum begins tc peak in angle at the energy source regions where 

Yj^(0)>O. 

This can be seen from numerical solutions to the steady state diffusion 
equation. In one such solution, shown in Figure la, the Langmuir energy has 
piled up in the forward propagating modes (e<n/2) which are in Cerenkov 
resonance with the positive slope of the beam. Figure lb displays clearly the 
angular gradient of the spectrum that may occur at some wavenumbers. As a 
result of this spectral behaviour, the weighted average defined by Equation (10) 
emphasizes the contribution from the forward angles, for which Yj^(e)>0, in 


relation to the wider angular range, for which Yj^( 0 )<O. By contrast, in 
the average, all angles are weighted equally. Therefore 




( 13 ) 


This inequality becomes stronger for more highly peaked spectra. In 
fact, <Y|^> may go through zero, in which case a {.hysical steady state 
solution (c£. Equation (12)) no longer exists, (i.e. the growing Langmuir waves 
cannot be saturated.) A discussion about marginal stability is deferred to 
the next section. 


Let us assume here that <Yjj>< 0. Then combining Equations (12) and (13) 
one has 




gfb),s(T)^S^M) 


(14) 


Where may be easily calculated from Eqs. (8)-(9) and evaluated for 
various parameters. turns out to be comparable in magnitude to M|^, 

In Table I values of Yf,^L y^^^, y^^^, and W^^"/T are reported 
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over a range of wavenumbers for the "nominal" beam and tail parameters. 

Uhen one compares the level of Langmuir waves determined by spontaneous 
emission without the beam to one notices that the angle-averaged 

spectral energy density is substantially enhanced due to the beam over a range 
of small wavenumbers only where the tail is not important) from Table I and 
Equation (14) it is seen that above kXjj = 5*10“^ the beam does not play any 


role> Thus, even if the bean>resonant spectral energy density, is 

enhanced by as much as 10^ through optimum choice of parameters, the effect 
on the total energv density is very weak. An upper bound to the enhancement 
of the energy density would be 

«J(kXp)(kXjj)V(2n^) • loVnxJ « 2‘10"^/nxj 

2 -2 

BO that the total energy density would be E /(^TinT) ® 5-10 /nx^ and 
the electric field only 1.5 uV/m for the solar wind parameters at 1 AU 
(n = 10 cm"^, Xjij = Bm, T = 13 eV). This is, of course, far too small to 
account for the observations. He are led to conclude that the rise of plasma 
waves measured when the beam arrives Is not due to an enhancement of the 
spontaneous Langmuir fluctuations as implied by Equation (14) but to growing 
waves which cannot be saturated by the processes considered. 

A comment should be made about the results obtained in a previous study 
(Goldman and Dubois, 1982). In that study the role of the electron tail was 
not considered, so that the beam enhanced modes could eKtend up to kx^ = 0.15, 
which resulted in an enhancement of the total energy density of two 
orders of magnitude, and an enhancement of one order for the electric field. 

3.3 Criterion for Steady State Solutions in the Presence of Isotropic Density 

Fluctuations 

When the density fluctuations are isotropic. Equation (6) is easily 
integrated, and the diffusion coefficient takes a form which is independent of 
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°lBO " (l/kXQ)2 1/k J dq C(q) 

" ”/12 «pe/(kXjj)^ q/k <6nVn^> 




Mhere q is a typical wavenumber associated with the fluctuations. Due to 
the k dependence, it is expected that unstable modes with a high phase 
velocity will be more quickly diffused out of the unstable region than will 
lower phase velocity modes. Thus, the quenching of the instability 
by the density fluctuations is expected to be much more effective for the high 
phase velocity waves. In fact, this difference in behaviour is even more 
pronounced because of the angular dependence of the Langmuir spectrum at 
different k. That feature should become clear for the reader from the 
following model. 


Let us consider the energy source and sink angular regions determined by 
the beam on a sphere of fixed radius k. For that purpose we rewrite Equation f8c) 
and Equation 9Cc) in terms of the variable x = cose with k considered as a 
parameter 


X exp [ -vJ/2ov 2 (x-Wpg/kv^)^ ] 


( 16 ) 


Sfb)(x) = 

X exp [ -vg/2av2 (x-Wpg/kVj^)^ ] 


( 17 ) 


The behaviour of the above expressions is characterized by three particular 



angles* ^ \s raaxinium, 

" '^pe'^^'^b ■*■ inantrnum, and Xj^ = “pe/*^''b “ 

whore is minimum. 

In the region Kq<x< 1 one has the power input 

P<„ = J W(x)dx + o’*" 

in "y 

’^o 

with o"^ = 2tj f S^^^Cx)dx>0 
«o 

while in the region -Kx<x^ one has the power output 

X 

^^out “ J fx) dx - o" 

with 0" = 2ti J^°S^^^(x) dx > 0. 

One necessary condition for stability is 

‘’in “ ‘’out • 


On the other hand, the flux of energy from one region to the other, 


» - 2rn-n2) D(.v) f = 2n D„ (i_,2, 


fl0) 


( 19 ) 


(20) 


( 21 ) 


must balance the input of power, 
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« 



( 22 ) 


which represents a second requirement for stability. 

« 

In order to form analytic criteria for steady state, we now model the 
spectral shape, according to our numerical experiments, by a Maxwellian 
centered on axis and a flat background (cf. the numerical result in Figure lb 
for an example)! 

W(x) = Wjjj + 6W exp [ - 5 ((x-D/CXq-I))^ ] , 

( 23 ) 


with the two free parameters adn tVi. Then we may rewrite Equation (20) as 
(the detailed calculations are left for the Appendix) 

SP < 1.17/ f 1-Bitp I-v?/26v2(1-x.)2] 1 _ 0.17 

L D 0 j j24) 

where the swelling factor SF is defined as the ratio of the peak to bottom 

values of the spectrum SP = 1 + aW/W . 

m 

Thus the very fact that the energy sink and source angular regions are 
finite on the sphere of radius k, and in a definite ratio, leads to a maximum 
bound of the swelling factor. As the marginal stability boundary is 
approached, the spontaneous emission becomes negligible and SF is equal to its 
maximum bound. Now, as the wavenumber k increases, x "pg/lfVjj decreases 
so that the maximum bound decreases too (cf. Equation (24)). The spectrum is not 
allowed to be as highly peaked for large k as it is for small k. However, a 
certain amount of peaking is required in order for the flux « to balance the 
input of power: Equation (22) leads to (cf. Appendix for details) 


r c^t? i o 


SP > 1 


0.4B D /w fvL,/-fliV)^x.fl+K-j . 

^ / C — Pe Dm, . Q ,P„.. -0.B5) 

^njj(Vb/6v) [ l-exp[-(v^^2av2)(l-Xjj)2] ] •> 


-1 


(25) 


Let UE define for a given wavenumber k, 


^raax ' ' (ti/Be)^/^ nbW^g/CkAv)^ 


(26) 


Combining the inequalities (24) and (25), an upper bound on for a 

ZliaK 

steady state solution to exist may be written as 


max ~ ^crit'’‘o' 


(27, 


with 


exp[-(v^/26v^)(l-x^)^] - 1 


0 


(28) 


From the viewpoint of numerical solutions of the diffusion equation it is 
possible, for a given k or x to increase yI^^ by increasing the 

O ninX 

strength of the beam n^ to a numerically determined maximum value, Pcrit**^o^’ 
where the marginally stationary boundary is reached and beyond which the 
spectrum grows without (linear) bound. Due to the exponential dependence in x^ 
displayed by the RHS of inequality (27) Pcj-itCx^) is expected to fall-off very 
rapidly with increasing k. 


In Figure 2 the exact obtained by the numerical solution of 

Equation (4) for given sets of parameters characterizing the beam has been 
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plotted as a function of kv^^/Upg s= l/x^. The agreement with the 

mt 

eKpression for given in (28) is satisfying. The rapid fall-off 

with increasing wavenumbers or decreasing phase velocities is evident. By 
comparing with the relatively slow k~^ decrease of one sees 

EQ&X ■ 

that, even if the small k modes with the supposedly largest growth rate are 

stabilized because, for those k, r„,.. for larger k the 

crit max ’ 

Inequality may be reversed so that the beam instability is then only partially 
stabilized. Two important points are worth stressing about those "residual" 
unstable waves. As compared to the fastest growing k-mode with phase velocity 
Wp/k “ Vjj-6v and e ° 0, they have a smaller growth rate given by Equation 
(26) anu are oblique, cose x^^ < 1 . 

Let us now turn to the role played by the nonthermal tail. In the first 
place, the associated damping increases quickly with increasing wave 
number (cf. Equation (Bb)). Due to that, the stability boundary. 
takes a V-shape as plotted in Figure (3). The second effect caused by the tail 
is less transparent. It may be understood however following the same line of 
reasoning as before. On a sphere of radius k the nonthermal tail reduces both 
the energy source region and the strength of the source t correlatively the 
energy sink tregion as well as the strength of the sink is increased. This 
gives rise to a greater peaking of the spectrum. The flux of energy from the 
source area to the sink area can then be more important. Therefore, the input 

of power can be raised, so that is upshifted. For the plot we have 

used a diffusion coefficient = a-lO (kXjj) j it corresponds to 
a fluctuation level of 8n/n = 5-10"^ in the range q"^=100 km, which is 
low compared to levels measured by ISEE. The parameters of the tail have been 

fitted from the non-Maxwellian component observed by Lin et al. (1981), 

n,p(Vj,,Vg)^ = 20. In that situation a beam of parameters nj^=5>10"^. 


Vjj/Vg*90, AV/Vg»20, is completely stabilised (cf. In dashed 
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line). Yet, were the beam 40t denser, n^ey.io , a small range of k 
would be unstable (dotted line). 

In general, three different regions of k-space must be considered, 
small k region, k<k^, Is stabilized by the diffusion, k^ may be 
expressed by using Equation (2B) for evaluating Vit the condition 

Let Xj^ = ’^pe'^^l'^b* weans of Equation (26) one may obtain the implicit 
for 

= 1 - /2av/Vb {Ln[l + 0.58 Do/Wpevg^v/v^ 1/n^^ Xid+Ki)] ) 


In the high k region, k>k 2 , y^^^ is stabilized by k 2 is defined 

by 


v^’'^(k2) 




or, using Equations (Bb) and (26) 

k2 Xjj = (0.10 n^jVg/njV^ 


In the 


equation 


(291 


(30) 


In the intermediate region, kj<|c<k 2 , there might remain weakly unstable modes 
so that the beam instability is only partially quenched and therefore slowly 
growing electric fields are present. Then, either the instability saturates 


because of some weak nonlinearity such as decay off ions, or simply because 
after a Mhile (t lo min) the waves enter in Cerenkov resonance with the 
negative slope of the beam due to the advectlon of the latter. Depending on 
the actual parameters, the relative position of the two curves and 

^max change so that various growth rates, and therefore levels 
of measured electric fields, are expected. In Figure 4 we display the change of 
^crit about by: a) an Increase in the level of density fluctuations; 

b) a downshift of the beam velocity and c) a decrease in the strength of the 
tall. For comparison the of Figure 3 is indicated with a dashed line. 

We are now ready to write a sufficient criterion for linear saturation by 
letting kj^>k 2 and using Equation (29) 

kjv,-/2kjiv {Ln[U0.58 < “pe 

“pe"b^'<2^DJ ,,,, 


where k^ may be computed via Equation (30). 

It is important to note that when the strong level of density 

fluctuations actually measured by ISEE is used to evaluate the coefficient 
-B 

(D^ o 10” above inequality is satisfied for any conceivable 

beam associated with Type III bursts. Thus, the instability is definitely 
saturated linearly by the diffusion process. 

3.4 Anisotropic Density Fluctuations; the Quest for Instability 

At the end of Section 3.2 it was shown that the level of electric field is 
hardly raised above background when dealing with steady state solutions to the 


diffusion equation. On the other hand, it has just been claimed that the 
density fluctuations measured by ISGE Here so strong that a steady state 
solution is justified. Therefore, He are faced with a dilemmai either the 3 
hour measurements on September B, 1976 by ISEE Here quite exceptional and not 
at all representative of the usual solar wind, or something is basically wrong 
in the model developed. Me suspect that the difficulty lies with one 
assumption which is not supported by observational evidences; the assumed 
isotropy of the density fluctuations. 

Let us assume that the density fluctuations split into two regions fcf. 
Figure 5a) 

C, , jcosil > a 
CCe) = " _ 

Cj^ I cos© I < a 0 < a = cos©^ < 1 

Then Equation (6) may be easily inteorated to yield 

fl(6) = H(a-sin©) Cj^ti/2 

+ H{Ein©-a) { Cj |Ti/2+(Cj^-Cj | ) [ Arcsin(a/sinQ)-a/sin0(l-a^/Bin^©)^''^ ] j 
In case of isotropic density fluctuations with the same energy one obtains 


®iso “ 

Therefore, the relative reduction of the diffusion coefficient 
may be evaluated at some angles. 


Let us imagine for example that the fluctuations are concentrated in a 
plane perpendicular to the beam axis. C|j = o, a << 1 (cf. Figure 5cK Then 
afTT/2)/0jg^ = 2 <x^/ 3 << 1, namely the forward propagating modes are 
nearly isolated from the backward ones. The forward modes are generally 

I 

destabilized by the beam whereas the backward ones are always dampedso that 
the unstable domain of phase space is cut off from the one of absolute 
damping, and the beam instability may not be expected to be completely 
quenched. On the other hand, if the fluctuations are concentrated along the 
beam axis. Cj^ = 0, a i 1 (cf. Figure 5b), one has Qfo)/ajgQ == 0. namely 
the parallel modes do not experience any diffusion at all and growth persists. 

Of course, from a physical point of view, since the beam axis coincides with the 
magnetic field axis and because the diffusion of the plasma particles is 
inhibited by the magnetic field, the inhomogeneities in density are more 
likely to be oriented across rather than parallel to the axis. 

In view of that we introduce the following model of anisotropic density 
fluctuations 


C(q,e) = p(q) a/u (cos^e+a^)“^ a << 1 


(32) 


where p(q) is determined from the ISEE 1-2 experiment (Celnikier et al., 19B3) 


p(q) = 2.5 10"® 9.10'® < qX^j < 7-10"^ . 


By substituting Equation (32) into Equation (£) one obtains 


Q(e) = 7.6 10"^° a/sin^e [l - a/(a^+sin^e)^/^ ] 


Thus the diffusion coefficient (Equation (5)) reads 


I 


D(k,©) = 6.4 10"^ w (1/kXr,)^ a/sin^e [l - a/(a^+sin^Q)^/^ ] . 

pa u I. (33) 

The Pokker-Planck equation haa been solved numerically using this diffusion 

coefficient. In Figure 6 it may be seen that the anisotropy in the density 

fluctuations leads to a concentration of the Langmuir energy in the forward 

propagating modes. Now, if the anisotropy is pronounced enough, the flux of 

energy at large angles between the source and sink regions may be too small to 

balance the power generated. The largest angle we have to consider is given 

by cose = Wp/k 2 Vjj where k 2 is defined in Equation (30). On the other hand, 

•*6 

it has been shown in Sec. 3.2 that with the parameters, = 10 , 

Vb = 90 Vg. av = 20 v^ and n^fv^/Vg)^ = 20 a diffusion coefficient 

D, B B-10"^^(kXn)"^ is too small to provide enough diffusion to 
iso u 

quench the instability. Therefore, we may expect to recover this situation 
for a < 8.10"^^/6.4 10“®- (l-(Up/k 2 Vj 3 )^) = 10‘^. The numerical 
integration of the Fokker-Planck equation indeed shows that for a < 5- 10"^ 
the fluctuations are unable to quickly remove the Langmuir energy from the beam 
unstable region of the phase space, and the wave develops unstably in spite of 
the strong general level of density fluctuations. Only in this case, may 
electric fields in the mV/m range be expected (consistently with observations). 
Thus, although measurements of both electric fields in the mV/m range and electron 
density fluctuations in the 10"^ range have not been obtained simultaneously, 
our results strongly suggest that such density inhomogeneities of the solar 
wind must be elongated along the magnetic field. 




4. Conclusion 


This study of the role played by the large scale density fluctuations of 
the solar wind in relation to Till electron streams leads to the following 
conclusions. 

1) It is possible to understand the rather low level, intermittent 
character of the beam-excited Langmuir waves observed in situ, by considering 
the background density fluctuations. In an environment where the latter are 
strong and isotropic, the growth of Langmuir waves is quenched and the 
electric field is hardly raised above its background value. However, when the 
density fluctuations are smaller, slowly growing Langmuir waves are present, 
which may result In an electric field in the mV/m range. 

2) Another explanation for the mV/m Langmuir fields observed is as 
follows 1 If the relatively small scale density inhoraogeneities of the solar 
wind (Sn/n " 10 , q " 10 kms) are assumed to be highly anisotropic 

and elongated along the magnetic field in a way consistent with the large 
inhomogeneities (8n/n ~ 1, q"^ - 10^ kmsl invoked for ducting the radiowaves 
through the corona, then we have shown that they are relatively ineffective in 
scattering (diffusing) unstable Langmuir waves into stable regions of k-space. 
In this case, one may turn to nonlinear and/or quasilinear processes for an 
explanation of the saturation mechanism. 
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Append! K i Evaluation of o, with a model of spectrum 

a) PluK of energy 

Substituting the ansatz (23) into Equation (21) one obtains 
® « 2tr//e Dq (VjjVgl^ *‘o 


b) Input/Output of power Pj^, p^^^ 

In order to simplify the integrals in Equations (10)-(19) we replace the 
Maxwellian shape of the spectrum by a box with the same area; 


W(x) s 


W_ 


m 


W. 


m 


+ O.B5 aW 


for x>Xg 
for K<Xg 


(A2) 


The integrals are then straightforward, 

Pin «= o'*" + ✓2n^/^ njj { 1-exp [ -vJ/2av^(l-Xg)^ ] j (Wjjj+0.B5aW) 

”out ' "b {l-«'<p[-Vb/2bv2(l+«^)2] ] 

Using these expressions, the condition (201 may be rewritten in the limit 
Vj^/av»l, Xg>0, 

(1+0. B5 aW/W^ { 1-exp [ -vJ/2av 2 d-x^)^] } 

«= 1 - (o‘*’+o") k^awjj/(^2Ti^'^^ “pe"b^’ 

Bence the inequality (24) where the equality is reached when the spontaneous 
emission terms become negligible. The other inequality (25) may be 



obtained in the same way by using Equations (Al) and (A3). 


Fig. 1 Langmuir spectrum when the diffusion Is sufficiently strong to quench 


the beam inbtabllity by diffusing t.ie waves out of the unstable 
regions of phase space into stable ones, yet not strong enough to 
prevent a pronounced peaking of the spectrum in the energy source area. 

a) Contours of (unit Tg) in the k-e plane > the Langmuir energy 
piles up in the forward propagating damped modes (0<n/2) which 
are in Cerenkov resonance with the positive slope of the bebm and 
diffuses into the backward damped propagating modeu (e)n/2'>. 

b) in unit of Tg at hVjj/«pg = 1.G2; note the substantial 
angular gradient for 

The parameters used are: n. = 5- 10"^, v. » go v„, ov »= 20 v . 
n^(V/j./Vg>^ = 20, and isotropic density fluctuations with 
5n/n = 5. L0~^, 

Fig. 2 Marginal stability boundary of the beam growth rate as a 

function of kVj^/Wpg = (cf. Equation (20)). If, at a given 
kVjj/Wpg, the growth rate of the Langmuir waves maximised over 
the angle 9, = 0.38 nj^/(k&v)^, is smaller than 

^crit instability at that wavenumber is quenched. 

In the plot shown, the density fluctuations are chosen isotropic 
with 5n/n = 5-10“'^ and the electron distribution function consists 
of a Maxwellian background without tail and of a beam with varying 
parameters (solid line: v^ = go Vg. tv = 20 v^} white circles: 

Vj^ » 90 Vg, ov «= 30 Vg! crosses: Vj^ = 60 Vg, 6v « 20 Vg). 




The rapid fall-off of is contrasted with the slow k"^ 

decrease of indicated in dashed line for a beam of 
max 

nj^ = 5-10 , Vjj B 90 Vg. 6v « 20 Vg. Thus, due to the diffusion, 

waves with the fastest growth rate may be stable while waves with bigger 
wavenumbers and slower growth rate are unstable. 

Fig. 3 Marginal stability boundary as in Figure 2 except that a 

nonthermal tail in the electron distribution is considered. The 

parameters of the tail have been fitted from the nonthermal component 

observed by Lin et al., 1981, n,j,(v,p/Vg)^ - 20. In that situation, 

a beam of parameters nj^ = 5-10"^, v^^/Vg = 90, Ov/Vg = 20, 

is completely stabilized (cf. 7^*^^ in dashed line). Yet, were 

loaK 

-7 

the beam 40t more dense, nj^ «= 7.10 , a small range of kx^ 

would be unstable (dotted line). 

Fic.^_i for various parameters. The marginal stability boundary 

^crit Figure 3 is indicated in the dashed-line for comparison 
purposes. 

a) Vj^ = 90 Vg, ov = 20 Vp. ny(v,p/Vg)^ = 20, 5n/n « 1.5 10"^ 

b) same as for a) except v^^ = GO Vg 

5 

c) same as for a) except n,pfv,j,/Vg) = 2. 

Pig. 5 A simple model of anisotropic density fluctuations. The function 
C(e) (cf. Equation (6)) is assumed to be at two levels. 

Fig. 6 Spectrum of density fluctuations a) and corresponding spectrum of 

7 7 

Langmuir waves b). On panel a) the equilines of Log (q (6n/n) 1^^) 
are plotted in the q-e plane. The density fluctuations are assumed 


to be strong and highly peaked in the plane perpendicular to the 
beam or magnetic field asis (c£. Equation (19) with a ^ 7-10"^). 
On panel h) the equilines of Langmuir energy are plotted in the 
k-e plane for the "nominal'' beam-tail parameters (nj^ » lo“®, 

Vjj a 90 Vg, Av a 20 Vg , n.j.fv.p/Vg1^ 20). Note the substantial 
concentration of energy in the forward propagating modes. 


Table for Detailed Balance 


"nominal" parameters 
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Abstrac_t 


We consider the equilibrium states described by a damped and driven 
kinetic equation governing the evolution of the spectrum of radiation in a 
stationary Langmuir-turbulent plasma. Both Langmuir and transverse spectra 
are assumed to be one-dimensional in wave-number space. A source of radiation 
at the plasma frequency and a uniform rate of dissipation at the higher 
harmonics are assumed. The radiation is characterized by an effective 
temperature, z, proportional to the Langmuir energy density and inversely 
proportional to the dissipation rate of the transverse waves. If z < 1 the 
equilibrium photon spectrum decreases with increasing frequency as a power law 
in z, If z > 1, photons scattered out of certain regions of phase space may 
return at a significnat rate, and the spectrum is found to have a global 
maKimum above the plasma frequency, The location of this maximum is 
proportional to y which may be many times the plasma frequency. 

r 

Applications to a laser-plasma experiment and to the solar-wind environment 
are discussed. 


I. Introduction and Hodal Kinetic Equation 


Although there a number of important applications, little is known about 

mechanisms for producing radiation at frequencies much higher than the maximum 

plasma frequency in non-magnetic plasmas. Radiation at high multiples of the 

plasma frequency has been observed in laser-pellet experiments,^ in 

relativistic beam-plasma interactions,^ and in low-energy beam-plasma 
2 A 

interactions. ’ It is also thought to play a role in astrophysical 

e 

phenomena.^ 

In relativistic beam-plasma systems, it has been suggested^ that the 
high energy of beam electrons could be tapped by Compton-conversion of a 
beam-electron and a plasma wave into a high-energy photon together with a 
lower-energy electron.'' Recent work shows that the growth-rate for this 
process is very low for a wide variety of Langmuir spectra.^ 

Radiation at twice the plasma frequency can be produced by coalescence of 

two Langmuir waves, and the process appears to be well-documented by examples 

g 

of radio-wave emissions from the solar corona.'^ Bowever, the emissivity is 
proportional to the square of the intensity of the Langmuir-wave spectrum, and 
the process is therefore not very efficient unless the Langmuir turbulence is 
quite intense (Langmuir energy-density on the order of the particle 
energy-density). Emission at the n-th harmonic of the plasma frequency by 
n-plasmon coalescence is proportional to the n-th power of the Langmuir wave 
intensity, so that this process would be still less efficient at weak levels 
of Langmuir turbulence. 
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Another elternative is to begin with radiation near the plasma frequency 
or its second harmonic, produced by instability or by external sources , and to 
allow it to Raman-scatter repeatedly from pre-existing Langmuir turbulence. 

At each step such scattering will produce frequency upshifts (anti-Stokes 
process) and downshifts (Stokes process — which is permissible as long as the 
shift does not go below the plasma frequency). The resulting radiation 
spectrum may exhibit frequency components significantly above the plasma 
frequency, 

This is the process we shall consider in this paper. For sufficiently 
weak photon damping rates, we find that multiple Raman scatter to the n-th 
harmonic is not proportional to the n-th harmonic of the Langmuir Intensity, 
so this process can be more efficient than the multiple plasraon coalescence 
mechanism under certain conditions. 

The multiple Raman scatter process works best in a hot plasma, with the 
radiation in contact with the turbulence for a sufficiently long time. This 
could be expedited in various ways: In astrophysical applications the spatial 

extent of the turbulence may be considerable. In laboratory plasmas there may 
be multiple reflections from chamber walls. 

Consider the scattering interactions of two species of waves illustrated 
in Figure 1. In a plasma, for example, the waves labeled by k could be 
transverse waves scattering off ion acoustic waves labeled by q. This process 
is commonly called Brillouln scattering. If, instead, the waves labeled by q 
are Langmuir waves then the scattering is termed Raman. In addition to 
higher-order wave-wave interactions not depicted in Figure 1, Brillouin and 


Raman icattering generally compete in determining the evolution of a 
transverse spectrum. (Have-particle interactions are ignored in this paper.) 
He assume the presence of an enhanced background of Langmuir turbulence, 
(rather than ion acoustic turbulenc)> ) , and focus our attention on the effect 
such turbulence has on the evolution of an externally produced electromagnetic 
(transverse) spectrum. 

He take the Langmuir waves to be statistically stationaryi If W^(k) 
denotes the ensemble-averaged energy density of Langmuir waves at wave-vector 
then H*'(]c) is assumed not to change in time. For example, the Langmuir 
waves may be maintained in steady state by a source such as an electron beam. 
Zn any case, we ignore the back-reaction of the transverse spectrum on the 
Langmuir spectrum. He are justified in doing so provided the transverse 
spectrum is far less energetic than the Langmuir spectrum responsible for the 
scattering. To summarize, we~are considering a plasma in which the 
ion-acoustic, transverse, and Langmuir wave energy densities are ordered as 
follows. 


H^® << h’' « H^ 

He also rely on the second inequality to justify our neglect of the 
stimulated scattering of transverse waves in Sections II and III. This 
process is described by terms in the kinetic equation that are proportional to 
(H^) , whereas the terms that we retain in the equation are proportional to 
H^H^. In general, if stimulated scattering must be included in the 
evolution of the transverse spectrum then so must the evolution of the 
Langmuir wave spectrum, since in this case H^ " H^. (He briefly describe 
equilibrium and temporal evolution in the presence of stimulated scattering in 
Secs. IV and V. ) 


The scattering that we consider here is resonant, that is. 


k » k’ 1 q 


(2a) 


and 


u(k) = w(k’ ) + Up 


(2b) 


are enforced in the kinetic equation (see Pig, 1). Here Wp denotes the plasma 
frequency and we neglect the thermal dispersion of Langmuir Haves, 

Our final assumption is that the Langmuir spectrum is confined to only 
one direction in g[-space, i.e., that it is maintained by a narrow source. 

With our source of radiation taken to be in this same direction, the kinetic 
equation describes transverse-wave evolution in a one-dimensional k-space and 
therefore ignores effects such as angular diffusion. 

Prom the dispersion relation for transverse waves 

u(k)^ “Up + c^k^ 

and using Eq. (2a) and Eq. (2b) it is easily seen that only Langmuir waves with 
wave numbers in the interval 






contribute to the scattering. We ignore any variation of W^(q) on this 
interval (see Pig. 2), and simply assume a flat spectrum. 


According to Equation (2b), only transverse waves whose frequencies 

differ by a multiple of the plasma frequency will interact with one another. 

Therefore with a source of radiation at the plasma frequency, the transverse 

spectrum (i.e., the ensemble-averaged energy density) will evolve only at the 

harmonics of u„. Let n be the harmonic number, 

P 

n E «(k)/Wp, 

Then the usual “Golden Rule” kinetic equation® describing this scattering nay 
be written in the following form. (We outline a derivation of this equation in 
Appendix A. ) 


^ c _w (A? bJ) U„/n 

dt u n n' 

+ W (aJ n/(n+l)2 + bJ n/(n-l)2 U^_^) 

+ V (U^/n) (aJ Un+i/(n+l)^ - bJ U„^^/(n-l)2) 

■ 2 ^n ^n + (3a) 

Bere, is the transverse-wave energy density at the n-th harmonic normalized 
to the ambient thermal energy density of the electrons, 4'nn^Tg (n^^ is the 
average plasma density, and is the electron temperature), and time is 
measured in units of the plasma frequency, Wp. The coefficient W is 
proportional to the Langmuir energy density, and is defined as 

W E L (v^/c)^ o W^ , 


(3b) 


fege- B 


where is the thermal velocity of the electrons i c is the speed of light in 
a vacuum; is the total energy density of Langmuir waves normalized to 
^^"o^e’ is > numerical reduction factor that selects, from a more 

general Langmuir spectrum, that portion of the spectrum which is contained in 
q-space between v^/c and /SUp/c and therefore active in the scattering 
considered here: 


« < 0.7. 


Equality holds if all of the Langmuir energy is contained between u^/c and /Su^/c 


(see Fig. 2). 


The scattering coefficients are defined as follows.^ (See Appendix A.) 


A^ E /[(n+1)^ - 1] - /(n^-l) 


(3c) 


n-1 


n > 2 
n = 1 


(3d) 


Equation (3d) expresses the result of detailed balancing. Using it one can 
easily show that in the absence of sources and dissipation (y^ >= = 0) 

the total energy density, ^ Ujj, is conserved. 


The coefficient of the terms which describe stimulated scattering is 




Therefore, roughly speaking, the ratio of the stimulated to the linear 
scattering rate (M^U) is 


PESe 9 


«'\'VU/W . 

When this ratio is much less than one, the stimulated scattering is ignorable, 
and the evolution of the spectrum is linear in 1). 

For simplicity, we assume a source term for radiation at the fundamental 

only, 


Throughout this paper we take the source at Wp only, but it is straightforward 
to generalize our results to include different source configurations. 

An independent source of radiation at 2Wp (n=2) is the coalescence of two 
Langmuir waves into a transverse wave. An important question is whether or not 
the coalescence process produces a stronger source of radiation at 2Up than 
does Raman scattering of a given source at Vp into that second harmonic. In 
fact, multiple Langmuir wave coalescence into radiation at the nth harmonic of 
the plasma frequency generally should be compared to the nth-ordet radiation 
produced by multiple Raman scattering. Bowever, in the one-dimensional k-space 
we are assuming the matrix elements for coalescence vanish so that this process 
is absent. (See Ref. 9 pg. 354, for example.) 


d = 


We take all harmonics to be dissipated at a constant rate, y 


Pege 10 


I 


n « 2,3, 


(4a) 


Lose o£ energy at n«p by convection out of a bounded plasma occurs at a rate 
equal to the group velocity at nWp divided by the size of the plasma, so the 
dimensionless radiation dissipation rate is taken to be: 


y 


d 



(4b1 


where R is the characteristic length of the plasma. Ue have approximated the 
group velocity by the speed of light, c (appropriate for the higher harmonics): 
and, again, all rates are measured in units of the plasma frequency. 

With this scenario of injection at the fundamental, Equations (3) predict 
a cascade of energy up to higher harmonics via Raman scattering. Our primary 
concern in this paper is with the equilibrium spectra that result from this 
cascade. Clearly such equilibria must exist if the damping rate at the 
fundamental, is positive (incoherent source of fundamental radiation)! 

In this case all harmonics are dissipated exponentially fast, whereas the 
injection of energy at the source proceeds only linearly with time. If the 
damping rate at the fundamental is assigned a negative value, to simulate a 
linear instability at Up, then the cascade either culminates in a stationary 
state or else the energy in the spectrum diverges with increasing time. In 
Sections II through IV we assume that a stationary state is eventually reached 
by the solution of Equation (3a) for the particular choice of y-^. In 
Section II we restrict our attention to equilibria for which 

«xW^ » U^, n = 1.2,... 





fio that the linearized kinetic equation (i.e., Eq. (3a) without the stimulated 
scattering terms proportional to V) adequately decribes the stationary state. 
Under this assumption it is possible to find analytic approximations to the 
equilibria using a continued-fraction representation. That analysis is carried 
out in Section III. In Section IV we present examples of equilibrium spectra, 
including a case in which the nonlinear terms (stimulated scattering) are not 
ignored. Unbounded behavior is considered in Section V. Concluding remarks 
appear in Section VI. 


II. Equilibrium Spectra 


Uith set equal to zero for all n in the linearized verison of 
Equation (3a), we are faced with the task of solving an inhomogeneous system 
of linear, algebraic equations. The number of these equations is, in 
principle, infinite, although we expect the spectrum to fall off rapidly at 
sufficiently large n to permit a truncation. Rather than impose an a priori 
truncation, we proceed, inifitead, to solve the equations recursively, writing 
each as a product of unending continued fractions. This will yield the 
exact solution for all n. With one very good approximation to the scattering 
coefficients, the continued fractions are summable in terms of Bessel 
functions (cf. Sec. 111). This method is applicable to other discrete 
systems enjoying linear nearest-neighbor interactions, although summability 
depends crucially on the form of the coefficients. 

For sufficiently weak photon dissipation rates, certain solutions we 
obtain clearly demonstrate scattering to high harmonics of with 
efficiencies well above unity (see Fig. 3). This is contrary to the popular 
belief that the spectrum should in any case fall off as The 

reasoning which leads to that conclusion ignores the scattering back down to 
lower frequencies described by the "A” coefficients in Eq. (3a). The 
scattering down results in an accumulation or “bottle-necking" of transverse 
energy at intermediate frequencies that is responsible for the observed 
extremum in the spectrum that we will approximate analytically in Sec. III. 


For convenience, we introduce 


and write the stationary, linearized version of Equation (3a) in the following 
form. (Me suppress the exponents on the scattering coefficients.) 

t WA„ + MB„] P„ - MA„ - HB„P„.i = S Bj „ 

i o J 

For illustrative purposes, we truncate the system to only the fundamental, 
and its first harmonic, P 2 . Than we can easily solve Equation (6)} 


+ 


“ (4tt^ 4- WA, HAy) 


W2 aJ 


and 


Pi 


4^2 4- MAi WAp 


UA, 


Retaining three modes, the solution isi 
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2 v^+WA^ + 


1 

(4 v 2 +MAi+MA 2^ 

MA2 


(GY3+WA2+WA3) 
1 . 




Pj . Fj 


4Y2+MA1+WA2 

WAl 


- ( 6 Y 3 +WA 2 +WA 31 

WA^/Aj 


and 


P3 = F2 ’ 


GY3+MA2+MA3 

MA^ 


il) 


As we increase the number of modes retained in .:he truncation, the following 
pattern emerges. 


Let i. 
t 


be the continued fraction, 





E 


i 
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ffia) 


Nith elements defined by 


and 




’m 

f 

(A, E 01 

“ ‘^^^e+j-i 3-1 

3>1 


MA^ 

£>0. 

3=1 

1 

£o0, 

3=1 . 


Then 


k+1 


s n. 


«=0 


fS) 


Equations (B and 9) are the exact solution of Eq.(6) and give us, using 
Eq. (5), the exact stationary states of the linearized kinetic equation. 

Uriting the solution in this form is particularly useful because of the slow 
variation of with n for n 2. 3- (Notice from Eq. (3c) that is 

approximatly equal to one if n 2 3.) Ue were careful to factor this slow 
dependence out of the scattering terms in the kinetic equation for the following 
reason. 

With our assumption that the are constant (« y^) for n 1 2 and 

the approximation of the A^'s by unity for all n, we can find simple 
closed-form expressions for all with n 2 2 in terms of P^. Therefore, we 
will have approximations to all c.f. eq. (5), and, hence, to the 

efficiency of the net up-scattering in the stationary state. Notice from Eqs. 


(8 and 9) that nc do not need to specify or S In order to determine this 
efficiency. These parameters serve only to determine the overall size or 
total energy density of the spectrum and do not effect its shape. Of course 
assume that the particular choice of Vj leads to a stationery spectrum, 
although for sufficiently large, negative values of the behavior is 
unbounded in time. See Section V. 
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III. Analytic Approxinations 

Kith all A^'s approximated by unity and, for the case of a constant 
dissipation rate at the higher harmonics. 


^n * 


nil 
n 1 2 , 


He find from Eqs. (6a thru c) that the continued fractions have the following 
simple form if £1 1. 




(€ > 1 ) 


Z i K/Vjjj 


This continued fraction is summable as a ratio of Bessel functions,' 




where denotes the Bessel function of the first kind of order 

It follows from Eq. (9) that 




(12) 
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•nd thus from Sq. (5) that 


(13) 


Notice the recursion relation implied by, and implying:, Cq. (10): 


= 2(€+z)/z - 1/f^^j 


(14) 


If we multiply this equation through by and use Eq. (9? Ms disjoyer that 


t+1 


+ F 


e-1 


2U±zI p 

z e 


This familiar three-term recursion relation is just Eq. (€) (n>l) and is 
satisfied by Pg(z) » among a great many other choices. 

Bowever, the continued fraction expansion, Eq. (10), specifies F^ exactly 
as Eq. (12). 


Using simple properties of the Bessel functions it is easy to see that 


0 < J 2 ^^(z)/J^^l(z) < 1 if z > 0. 

Therefore, Eq. (13) makes sense (i.e., U^ xs positive for all n) provided 
only that is positive. Requiring to be positive, in the presence of 
a linear instability at the fundamental (y^ < 0) enables us to locate the 
threshold for saturating the instability using o'’r exact expression for F^^, 
Eqs. (B and 9). We have observed thot negative *;5lues of imply that the 
numerical solution of the kinetic equation is not bounded in time. We discuss 
unbounded behavior in Sec. V. 
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Me win •ssune that is positive and conrider some physically 
interesting limiting forms of our expression for the equilibrium radiation 
spectrum. The accuracy of our formula, Eq. (13), will be tested by comparing 
its predictions with numerical solutions of Eqs. (3). 


IV. Examples 


It Is noteworthy that larger values of z result In more efficient 
up-scattering. Mhen the scattering rate exceeds the dissipation rate at the lower 
harmonics (z > 1) significant up-scattering occurs, producing an extremum in the 
spectrum at a harmonic of the fundamental. Uhen z < 1 no such extremum 
occurs, and most of the total energy resides at the fundamental. 

Me now consider the two limits: 

A. z >> 1: Strong Un-Scattering 

Sere we assume that most of the energy in the spectrum resides at n << z. 

Then using Debye's asymptotic formula for the Bessel functions^^ and retaining 
only leading terms in the small parameter, n/z, we obtain the following 
expression for the equilibrium spectrum in this case. 

U A exp (1-n^''^^] 

n i i z 

In Fig. 3 we plot the equilibrium spectrum obtained by solving the 
linearized Eqs. (3) numerically, in the limit of long time, with z = 10^ and 
Notice that the energy densities are small enough to justify our 
neglect of the stimulated scattering. (Here, as in Cases B and C below, we 
take the dissipation to be due to convection, e c/RWp, and ignore 
collisional damping, for example.) 


Agreement between the numerical solutions and the predictions of our 
formula, Eq. (15), is quite good and Improves at the higher harmonics. In 
particular, we call attention to the maximum in the spectrum at n » 12. The 
formula predicts a maximum at 


n 


m 


f JL.)2/3 -1/3 


which is equal to 11.5 if z is 10^. 


Certain of the parameters used in Pig. 3 are characteristic of the solar 
corona, although it is important to note that Langmuir wave energy densities 
as high as those we are assuming («W^ ^ o.7) have never been observed in 
the solar wind, where oW^ <. 10"^ is the rule. Such low levels of 
background Langmuir turbulence will not produce the strong up-scattering 
illustrated in Pig. 3, in a plasma the size of the solar corona. 

It may, however be worth noting that larger (hypothetical) stellar coronas 
could demonstrate strong up-scattering of fundamental radiation, even at weak 
levels of Langmuir energy density. To produce Fig. (3) with aW^ *» 7x10"® 
we would require a corona 10^^ metres in radius, or 10^ times the size of 
the solar corona. (See the caption to Pig. (3), and recall that z " 
where R is the spatial extent of the plasma.) 

Similar conclusions have been reached in a study of the diffusive limit 
(i.e., the high frequency limit) of the kinetic equation by Colgate and his 
co-workers. As noted by them, the discovery of a maximum in the radiation 
spectrum at a harmonic of the plasma frequency may have important consequences 
for the theory of radio emission from qunsi-stellar sources. 


B. 7 << 1: Laser-Plasma Interaction 


When 7 << 1, we may ignore it in the order of the Bessel functions. 
Expanding about 7 «= 0 and keeping only leading terms in z we obtain the 

following power-law formula for the equilibrium spectrum in this case. 


U A 


Ujn“ 


• 2 ^ 


n-1 


JL 

n! 


(16) 


Por a physical example of this limiting form of the spectrum we consider the 
observations of back-scattered light in a laser-target interaction reported by 
Burnett and his coworkers. In this case a strong source propagates through a 
plasma of small extent. A typical photon remains in the plasma for a relatively 
Short time and so is not likely to be up-scattered. The result is a spectrum 
that is strongly peaked at the plasma frequency. A value of z equal to 10"^ 
is appropriate for their experiment, and we have taken 

In Fig. (4a) we plot the equilibrium spectrum obtained by solving the 
lineariced kinetic equation numerically. In Pig. (4b) we plot the 
corresponding result for the full kinetic equation, including tha stimulated 
scattering. Crosses indicate experimental results reported by Burnett et al. 

Notice that efficiencies oberved in the laser experiment are greater than 
those predicted by either the linear or nonlinear kinetic equations at high 
frequencies. The scattering at these high frequencies is off Langmuir waves 
with wave vectors clustered about Wp/c. Thus, enhanced levels of plasma 
turbulence near this wave number (i.e., that wave number excited by the laser) 
would account for the observed enhanced efficiencies, in the context of Raman 
up-scattering. 


The nonlinear theory predicts snaller efficiencies than the linear 
theory. The effect of stimulated scattering is to suppress up-scattering and 
to concentrate a greater fraction of the total energy at the fundamental, 

«p, compared to the predictions of the linear theory. This result may be 
inferred by comparing the signs of the two nonlinear terms in the full kinetic 
equation (3a) > radiation at n+1 causes to increase with times radiation 
at n-1 causes to decrease with time.' Thus, stimulated scattering 
transfers energy to low frequencies, and so supresses up-scattering. 

Our numerical results are compared with the prediction of our formula for 
this case, Eq. (16), in Table I. Agreement with the linearized kinetic equation 
Is excellent at the third and higher harmonics, as it should be, since our 
approximation to the scattering coefficients is very accurate at these harmonics. 

In the experiment, the laser-source is so energetic that the stimulated 
scattering cannot be ignored, and this implies that the back-reaction of the 
Langmuir spectrum to the transverse waves ought to be considered as well, as 
discussed earlier. However, the stationary spectrum, obtained by integrating the 
nonlinear kinetic equation (3a) numerically, is reached very quickly. Starting 
from no radiation initially, we find that a stationary spectrum is reached 
within five pico-seconds, which is much shorter than the duration of the laser 
pulse (two nano-seconds). This means that the predicted equilibrium Raman 
spectrum is observable during the experiment. But it also suggests that the 
back-reaction on the Langmuir spectrum may be ignerable, at least for tiroes 
of interest in this particular example. 


C. Coherent radiation frow relativlEtic electron beams 


Recently, radiation has been detected at high harmonics of the plasma 
frequency during intense, relativistic beam-plasma laboratory experiments 
The reported spectrum is relatively flat, continuous (i.e,, no spikes), and 

extends from « out to at least 7 w_. 

P P 

The lack of spikes in the spectrum may result from a broadband instability 
and/or source of radiation excited by the electron beam, having a range in 
fre»iuencles on the order of The theory could be genralized easily to 
include such a mechanism for Injection near Me expert the beam to excite 
a linear instability at the fundamental, for example, via the parametric decay 
of a Langmuir wave into an ion-acoustic and a transverse wave. 

However, to explain a flat spectrum out to at least 7 in the context 
of the present theory, we require a value of z at least on the order of unity. 

2 B W/y^ >. 1. 

The Langmuir energy density, is not measured in the experiment, so the 
parameter W is undetermined. To ensure the validity of the weak turublence 
theory used in this paper, must be less than one. Therefore W is bounded 
above as follows (cf. Eq. (3b)) 


where a is less than 0.7, as explained in Sec. I. 


If the dominant dissipation mechanism is tssumed to be absorption at, or 
transmission through, the chamber walls, then = c/WpR, where R is the 
effective length of the cha^nber. In this case we estimate the following upper 
bound for z. 


Typical measured values of <v^/c)^ and «p are 10“^ and 10^® s"^, 
respectively. It follows that R must be at least 4x10^ m to give us a value of 
z that is greater than one and so to admit the present theory as an explanation 
of the observations. The value of R depends on the rate at which radiation is 
reflected by the chamber walls. If the walls do not reflect the radiation at 
all, then R is the actual chamber length (1 m in the experiment). In the 
limit of perfectly reflecting walls, R diverges the plasma is 

infinitely long, in effect, and efficient up-scattering is possible. In fact 
is bounded away from zero by other dissipation mechanisms, such as 
collisional losses. 

It is unlikely that the effective length of the chamber is as large as 10^ 
times the actual chamber length, so it is doubtful that Raman scattering is 
responsible for producing the observed high«frequency radiation. It would be 
interesting, however, to compare experimental results for a variety of chamber 
wall coatings differing in reflectivity. 


V. Unbounded Behavior 


Here we briefly consider some interesting properties of unbounded 
behavior resulting from linear instability at the fundamental. Nith the 
addition of a growth rate at the fundamental < 0) it is possible for the 
energy density of the radiation to diverge in time, eventually reaching the 
nonlinear regime where the stimulated scattering plays an important role. As 
mentioned in Sec. IV, stimulated scattering tends to push the transverse 
spectrum back down to low frequencies and thus to decrease the net rate of 
dissipation. As the energy density diverges (Fig. (5a)) it does so faster 
than exponentially (Fig. (5b)) with an increasing fraction of the total 
energy at the fundamental (Fig. (5c)), 

Bowever, it is important to emphasize that even though the long-time 
behavior is unbounded, the spectrum may be acceptably stationary for some 
finite time of interest, depending on the particular problem. For example, 
if the source is pulsed (along with the growth rate) it may shut off before 
the stimulated scattering is activated. 

There is a very interesting short-time behavior that is generally observed 
in these “unbounded" examples which dramatically illustrates the effect of the 
stimulated scattering on the evolution of the spectrum. In Pig. (5d) we see the 
linear evolution responsible for up-scattering the radiation to a mean frequency 
of just under 5Wp after about 0.7 seconds. (Except for the growth rate at 
the fundamental, parameters are similar to those in the corona example of Fig. 
(3).) At this point, the spectrum is as shown in Fig. (5e) and is on the 


t 
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brink of activatin? the atimulated scattering. Uhen it does so, the nonlinear 
terms push the radiation back down into the fundamental so that the average 
freqeuncy swings back down toward u^. 

A growth rate at the fundamental does not necessarily produce unbounded 
behavior. The growth rate must be larger than a certain threshold value above 
which the linear instability cannot be saturated. The case depicted in Figs. 

(5) lies just above threshold. Of course the equilibria reached in spite of 
linear instability may lie well within the domain of the linearized kinetic 
equation. In this case, the linear analysis correctly predicts the 
characteristics of the spectrum. Me emphasize that the shape of the spectrum is 
independent of the source and growth rate at the fundamental (see Eq. 12, for 
example), these serving only to determine the overall size (total energy 
density) of the spectrum. 


IV. Concluding Remarks 


Me have obtained the stationary states of a linear kinetic equation that 
describes the scattering of a source of radiation by a Langmuir-turbulent plasma 
of finite spatial extent. Because the scattering proceeds via nearest-neighbor 
interactions only, it is straightforward to write the solution in terms of 
unending continued fractions. These fractions provide an efficient algorithm 
for calculating the equilibria as well as accurate approximations to those 
equilibria in terms of simple analytic functions. The generalization of this 
method to include the effects of nonlinear nearest-neighbor interactions 
(stimulated scattering) is currently under investigation. 
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Table 1 


al Formula b) Kinetic Equation 


1 

1.00 

0.45 

2 

1.43 

1.33 

3 

1.65 

1.61 

4 

1.81 

1.78 

5 

1.92 

1.91 

6 

2.01 

2.00 

7 

2.09 

2.08 

6 

2.15 

2.15 

9 

2.21 

2.21 

10 

2.26 

2.26 

11 

2.30 

2.30 


Table Ii 

“log^Q[Ujj^^/Un] approximation to the spectrum given 

by eq. (16), and b) the long-time asymptotic solution of linearized Eq. f3a). 
z » 10"^ and physical parameters are given in the caption to Pig, 4. The 
absolute error in both columns is ± 0.005. 
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Figure Captions: 

1) Diagrams ataman scattering, k denotes photonsi q denotes Langmuir waves. 

2) The (one-dimensional) Langmuir wave spectrum that we assume. ® (energy 
density of Langmuir waves at q)/(4nnpTg). 

3) vs n for parameters appropriate for a hypothetical stellar corona with 

a higher level of Langmuir turbulence than is observed for the solar corona. 
Tg « 10^ eV, ■ 6KlO^/sec, ^nn^Tg ■ 2x10’^ erg/cro^, R = 10^ m, 
z « 10^. Tjj « 5x10“®. S = 100. » 0.7. 

4) vs n for the parameters of the laser-target interaction 
experiment of Burnett et al.^^ Crosses indicate their results. » 10^ 
eV, V B 1.7BxlO^^/sec, 4Tm.T- = 2x10® Joules/cm^. z = 10"^, 

P U V 

R « 3x10“^ cm, « 5.6x10“^, S »= 638, * 0.7, 

a) Linearized kinetic equation prediction. 

b) Full kinetic equation prediction, including stimulated scattering. 

5) Unbounded behavior resulting from linear instability at si, , 

a) Total, dimensionless energy density (J] vs. time. 

b) The logarithmic derivative of the total number density, i n, vs. 
time, i.e., an effective growth rate. 

c) The fraction of the total energy density that is at vs. time. 

d) The average frequency, ^ 

e) The spectrum at the time at which the up-scattering is greatest. 


SB 


2>!lO"^ erg/cm^, 


Tg » 10^ eV, ttp - 
R = 10® n, z e 10^, Yj = Yj o 5K10"®. 

«M^ « 0.7. 


S “ 10, 


F'Epe 
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Appendix Aj Reduction of the Golden-Rule Kinetic Equation to One-Dimension 


In terms of number densities for the transverse and longitudinal spectra, 


and 


BthuJ. V " "k 


' q' = 

8771^-.' V " q 


the well-known kinetic equation^ for the scattering processes described in 
Fig. 1 is 


(277} ■ 


{I 


© ® © 


• 2 L L 

M~| 6(U) -w. -w ) [-n n +n, n -n n, ) + 
' k k-q p k q k-q q k k-q 


0 0 © 

^2 L L 

+ |m I 6(u) -oj.^ -ho ) [-n, n +n. ^ n +n, n. 1 
' ' ' k k+q p k q k+q q k k+q 


fAl) 


where 


M 


± 1 2 _ 2 . e .2 

= ” <s-> vcr 

e k k±q 


Here we neglect the thermal dispersion of Langmuir waves in using w alone for 

P 

their frequency. The polarization factors usually present in the matrix elements, 

have been taken to be a constant and are suppressed. That is, we assume 

an unpolarized sci-rce. In Eq. (Al) we are also suppressing the vector notation 

on the wave vectors Jc and £. Terms 0 through 0 give us the linearized kinetic 

equation? tetins (§) and (|) provide the stimulated scattering. (Notice that we have 

defined Fourier transforms in an arbitrarily large cube in x-space of volume V 
j T I 2 I Xi I 2 

so that E, and E have the dimensions of energy •volime. ) 
k q 




We assume that both the Langmuir and transverse spectra are strongly 
peaked £ibout a single direction in k-space and calculate narrow-angle 
averages of the number densities about that direc 'i. The result is expressed 
in terms of the transverse wave number density as a fu.iction of frequency, 


N^k) 


^ 1 ^ 


For example, consider term in Eq. (Al), Integrating over we 


f da f d^q ^ n^ f da 6 (w. -i 
J k'^ j 0.' q J k k 


W, , -Hi) ) 

k+q p 


To perform the integration over define 

L - K+a ' 

and assume that the only contribution to the integral comes from a very narrow 

neighborhood of 6^ = 0/ see Fig. 6. Then, ignoring the variation |^| that 

results from a variation in 6j^, we have 

^dCcosBf.) = kd{cos9 ) 
t K 

2 

from = q +k,*a» So the integral over is approximated by 




and we have 


j li-il" 5^ "|.j, m5)f (v“5%> • 

Here ri(?) = | and is a function of 1^| only; d^^ «= d^q, now, with k 

fixed. Integrating over the 6-fxuiction in (d^, and using 
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W- 

6((J “W -MU ) - 6(C-^*) 

2 2 2 2 

where c (^*) = (u, +o) ) - w , leaves us with 
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j®d1, -V nCf) ^ j dn^. |£.-k|^ n^._^ . 


From 


(h*C*)cos0^^ - “ k cos6^^_j^ , 

and neglecting d(I_C*~k|), we have 


li*-k 


d{cos9^J = 


d(cos0^^_j^) 


or 


ii*-k 




Thus 




TKC*) (t:*)' 


^*^?*-k "e*-k 


c'kto^C* ^ - 

K 

The contribution to N(tJjj) is obtained after multiplying this expression by 

k^/Vg(k) , 


where V (k) = du), /dk. 
g k 

Using 


2 

~r = ^ > 

g k g 


we have finally 
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as the contribution ths ; te.Tii ^makes to N(u)j^). 

Using the same tricks to evaluate the remaining five terms and assuming 

that 

W S f dn n^ 

J n U 


is independent of over 


Ann 


0) D 

— ^ < |u I < /3 — ^ 
c — — c 


we arrive at the following kinetic equation: 

A . . N(ti) ) 

N(uj ) + ^ wtik-rr + k-rrH(a)^-2o) )] 

c k 
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where H is the step function, 


|k-5‘l =i ItoV)* - t(V“p>"-"p'*l • 


and 
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N(tOj^) has the dimensions of (time) / (length) ' 


Introducing 


n = 03. /u , 
h P 


a dimensionless energy density 


Q 




N{w ) 

Utn) = . 

o e 
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and measuring time in units of OJ , the kinetic equation is written in the 

P 


following dimensionless form: 


U(n) = -W[A^+HB^] U(n)/n + 


+ W[A^ U(n+1) + HB^ U{n-D) 

(n+1) (n-1) 


^ V [A^ liiillil - HB^ , 

" (n+1)^ (n-1) 


where 
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Recapitulating 


U(n) = -T — ) = *:; — 
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Dissipation will add a term to the right-hand-side of Eq. (A2) 


-2Y(n)*U(n) , 




where -2y{n) is the rate at yhich transverse energy is removed from the system, 
normalized to the electron plasma frequency. A source of radiation at k will add 
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Fig. 5b 
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Fig. 5d 






